Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions.

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Acknowledging more biodiversity without more species.

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

Kukri snakes Oligodon Fitzinger, 1826 of the Western Palearctic with the resurrection of Contia transcaspica Nikolsky, 1902 (Reptilia, Squamata, Colubridae).

The kukri snakes of the genus Oligodon Fitzinger, 1826 reach the westernmost limits of their distribution in Middle and Southwest Asia (Afghanistan, Iran, and Turkmenistan), and the Palearctic portions of Pakistan. In this article, we review the systematics and distribution of the two species native to this region, Oligodon arnensis (Shaw, 1802) and Oligodon taeniolatus (Jerdon, 1853) based on an integrative approach combining morphological, molecular, and species distribution modeling (SDM) data. Phylogenetic analyses recover O. taeniolatus populations from Iran and Turkmenistan in a clade with the O. arnensis species complex, rendering the former species paraphyletic relative to O. taeniolatus sensu stricto on the Indian subcontinent. To correct this, we resurrect the name Contia transcaspica Nikolsky, 1902 from the synonymy of O. taeniolatus and assign it to populations in Middle-Southwest Asia. So far, Oligodon transcaspicus comb. et stat. nov. is known only from the Köpet-Dag Mountain Range of northeast Iran and southern Turkmenistan, but SDM mapping suggests it may have a wider range. Genetic samples of O. "arnensis" from northern Pakistan are nested in a clade sister to the recently described Oligodon churahensis Mirza, Bhardwaj & Patel, 2021, and are phylogenetically separate from O. arnensis sensu stricto in south India and Sri Lanka. Based on morphological similarity, the Afghanistan and Pakistan populations are assigned to Oligodon russelius (Daudin, 1803) and we synonymize O. churahensis with this species. Our investigation leads us to remove O. taeniolatus from the snake fauna of Afghanistan, Iran, and Turkmenistan, with the consequence that only Oligodon transcaspicus comb. et stat. nov. and O. russelius are present in these countries. Additional studies are needed to resolve the taxonomy of the O. taeniolatus and O. arnensis species complexes on the Indian subcontinent, and an updated key for both groups is provided.

Global Protected Areas as refuges for amphibians and reptiles under climate change.

Protected Areas (PAs) are the cornerstone of biodiversity conservation. Here, we collated distributional data for >14,000 (~70% of) species of amphibians and reptiles (herpetofauna) to perform a global assessment of the conservation effectiveness of PAs using species distribution models. Our analyses reveal that >91% of herpetofauna species are currently distributed in PAs, and that this proportion will remain unaltered under future climate change. Indeed, loss of species' distributional ranges will be lower inside PAs than outside them. Therefore, the proportion of effectively protected species is predicted to increase. However, over 7.8% of species currently occur outside PAs, and large spatial conservation gaps remain, mainly across tropical and subtropical moist broadleaf forests, and across non-high-income countries. We also predict that more than 300 amphibian and 500 reptile species may go extinct under climate change over the course of the ongoing century. Our study highlights the importance of PAs in providing herpetofauna with refuge from climate change, and suggests ways to optimize PAs to better conserve biodiversity worldwide.

Genome-wide single nucleotide polymorphisms reveal recurrent waves of speciation in niche-pockets, in Europe's most venomous snake.

Within the Balkan Peninsula, topographic and climatic agents have promoted biodiversity and shaped the speciation history of many ectotherms. Here, we targeted an iconic European reptile, the nose-horned viper species-complex (Vipera ammodytes), and explored its spatial and temporal evolution. We (i) utilized genome-wide single nucleotide polymorphisms to infer genetic structure and build a time-calibrated species-tree, and (ii) applied species distribution modelling with niche-divergence tests among major phylogenomic clades. Geographically structured genetic diversity was found. Cycles of recurrent isolation and expansion during glacial-interglacial periods led to allopatric speciation and to secondary contacts and formation of multiple hybrid zones throughout the Balkan Peninsula. Deep divergence is still detected among populations separated by old and imminent biogeographical barriers (Pindos Mountain Range, the Cyclades islands, etc.), but in most cases speciation is incomplete. At the other end of the speciation continuum, we recognize two well-differentiated lineages, currently lacking any evidence of gene flow; one is distributed in the Northwestern Balkans and the other in the Southeastern Balkans, further expanding into Asia. Despite their split 5 million years ago, there is no evidence of ecological divergence, as speciation probably occurred in niche-pockets of analogous environments. These two lineages probably represent different species, while V. transcaucasiana does not merit species status. By comparing the genomic phylogenies to an updated mitochondrial one, we propose an evolutionary scenario that resolves all mitonuclear conflicts, according to which the history of the V. ammodytes species-complex was shaped by complex processes, including a major event of introgressive hybridization with asymmetric mitochondrial capture.

A new, rare, small-ranged, and endangered mountain snake of the genus Elaphe from the Southern Levant.

The genus Elaphe Fitzinger, 1833 includes 17 species of charismatic, large-sized, non-venomous, Eurasian snakes. In the Western Palearctic, the genus is represented by three species from the Elaphe quatuorlineata group ranging from the Apennine peninsula to Central Asia. The southernmost population of this group is distributed in the mountains of the Southern Levant, with more than 400 km gap to other Elaphe populations. This population has been known to science for only 50 years and is virtually unstudied due to its extreme rarity. We studied these snakes' morphological and genetic variation from the three countries where they are known to occur, i.e., Israel (Hermon, the Israeli-controlled Golan Heights), Lebanon, and Syria. We used nine mitochondrial and nuclear genes, complete mitogenome sequences, and a comprehensive morphological examination including published data, our own field observations, and museum specimens, to study its relationship to other species in the group. The three currently recognized species of the group (E. quatuorlineata, E. sauromates, E. urartica), and the Levant population, form four deeply divergent, strongly supported clades. Three of these clades correspond to the abovementioned species while the Southern Levant clade, which is genetically and morphologically distinct from all named congeners, is described here as a new species, Elaphe druzei sp. nov. The basal divergence of this group is estimated to be the Late Miocene with subsequent radiation from 5.1 to 3.9 Mya. The revealed biogeography of the E. quatuorlineata group supports the importance of the Levant as a major center of endemism and diversity of biota in Eurasia. The new species is large-sized and is one of the rarest snakes in the Western Palearctic. Because of its small mountain distribution range, in an area affected by land use and climate change, the new Elaphe urgently needs strict protection. Despite political issues, we hope this will be based on the cooperation of all countries where the new species occurs.

Phylogenomics of Anguis and Pseudopus (Squamata, Anguidae) indicates Balkan-Apennine mitochondrial capture associated with the Messinian event.

A dated phylogenetic hypothesis on the evolutionary history of the extant taxa of the Western Palearctic lizards Anguis and Pseudopus is revised using genome-wide nuclear DNA and mitogenomes. We found overall concordance between nuclear and mitochondrial DNA phylogenies, with one significant exception - the Apennine A. veronensis. In mitochondrial DNA, this species forms a common clade with the earliest diverging lineage, the southern Balkan endemic A. cephallonica, while it clusters together with A. fragilis in nuclear DNA. The nuclear phylogeny conforms to the morphology, which is relatively similar between A. veronensis and A. fragilis. The most plausible explanation for the mitonuclear discordance is ancient mitochondrial capture from the Balkan ancestor of A. cephallonica to the Apennine population of the A. fragilis-veronensis ancestor. We hypothesize that this capture occurred only in a geographically restricted population. The dating of this presumed mitochondrial introgression and capture coincides with the Messinian event, when the Balkan and Apennine Peninsulas were presumably largely connected. The dated nuclear phylogenomic reconstruction estimated the divergence of A. cephallonica around 12 Mya, while the sister clade representing the A. fragilis species complex consisting of the sister species A. fragilis-A. veronensis and A. colchica-A. graeca further diversified around 7 Mya. The depth of nuclear divergence among the evolutionary lineages of Pseudopus (0.5-1.2 Mya) supports their subspecies status.

The genus Blanus (Reptilia: Blanidae) in Iraq: distribution, species identification and genetic affiliation.

N/A.

Genotyping the phenotypic diversity in Aegean Natrixnatrixmoreotica (Bedriaga, 1882) (Reptilia, Serpentes, Natricidae).

We examined the mitochondrial identity of Aegean Natrixnatrixmoreotica representing different morphotypes, with a focus on new material from Milos and Skyros. We found no correlation between distinct morphotypes and mitochondrial identity. Our results support that grass snake populations are polyphenetic and that southern subspecies, including island populations, show a higher variability than northern ones.

Appearances often deceive in racerunners: integrative approach reveals two new species of Eremias (Squamata: Lacertidae) from Pakistan.

Based on newly provided morphological and previously published genetic data, we describe two new distinctive sympatric lacertid lizards of the genus Eremias (subgenus Aspidorhinus) from the arid mountains of northwestern Balochistan Province, Pakistan. The new species, Eremias killasaifullahi sp. nov. and Eremias rafiqi sp. nov. are distinguished from all other species of the subgenus Aspidorhinus (E. afghanistanica, E. fahimii, E. isfahanica, E. kopetdaghica, E. lalezharica, E. montana, E. nikolskii, E. papenfussi, E. persica, E. regeli, E. roborowskii, E. strauchi, E. suphani, and E. velox) by unique morphological characters and genetic differentiation. In the molecular phylogeny of Aspidorhinus, both new taxa are well-supported lineages differentiated from other species of this subgenus by uncorrected p distances from 8.5% to 21.6%, respectively. Both new species belong to E. persica complex where E. rafiqi sp. nov. is partly similar in dorsal color pattern to E. persica but can be distinguished from this species by unique meristic and morphometric characters. Eremias rafiqi sp. nov. is found in Pakistan, Afghanistan, and eastern Iran. Eremias killasaifullahi sp. nov. can be differentiated from E. persica by its distant distribution, dorsal color and pattern, smaller size, and less number of gulars and ventral scales. Eremias killasaifullahi sp. nov. is only known from the type locality and represents the local microendemism, along with other endemic species of reptiles reported from this part of Pakistan. We, however, expect that E. killasaifullahi sp. nov. could have a broader range in northwestern Pakistan and southeastern Afghanistan, which should be an object of following investigations. Our data show that remote areas between Hindu Kush Mountains and Indus River need attention as they most probably represent possible sources of genetic and species diversity in the region.

Temperature regulation in the Balkan spadefoot (Pelobates balcanicus Karaman, 1928) at the beginning of nocturnal activity.

On land, the amphibians interact with the environment in a complex way-even small changes in the physiological conditions may significantly impact the behaviour and vice versa. In ectothermic tetrapods, the transition from inactive to active phase may be related to important changes in their thermal status. We studied the thermal ecology of adult Balkan spadefoots (Pelobates balcanicus Karaman, 1928) in northeastern Bulgaria. These toads spend the daytime buried between 10 and 15 cm in sandy substrates, and emerge after sunset. On the substrate, their thermal energy exchange is defined by the absence of heat flow from the sun. Secondary heat sources, like stored heat and infrared radiation from the soil play an important role for the thermal balance of the active spadefoot toads. At the beginning of their daily activity, we measured substrate temperature (at a depth of 11-12 cm), toad's surface body temperature, and also provided thermal profiles of the animals and the substrate surface in their microhabitats. In animals which recently emerged from the substrate, the temperature was comparatively higher and was closer to that of the subsoil on the spot. After that, body temperature decreased rapidly and continued to change slowly, in correlation with air temperature. We detected a temperature gradient on the dorsal surface of the toads. On the basis of our measurements and additional data, we discuss the eventual role of air humidity and the effects of surface and skin water evaporation on the water balance and activity of the investigated toads.

The evolutionary history of an accidental model organism, the leopard gecko Eublepharis macularius (Squamata: Eublepharidae).

The leopard gecko, Eublepharis macularius, is a widely used model organism in laboratory and experimental studies. The high phenotypic diversity in the pet trade, the fact that the provenance of different breeding lines is unknown, and that distinct Eublepharis species are known to hybridize, implies that the continued use of E. macularius as a model requires clarity on the origin of the lineages in the pet trade. We combine multi-locus sequence data and the first range-wide sampling of the genus Eublepharis to reconstruct the evolutionary history of the Eublepharidae and Eublepharis, with an updated time-tree for the Eublepharidae. Our sampling includes five of the six recognized species and additional nominal taxa of uncertain status comprising 43 samples from 34 localities plus 48 pet-trade samples. The Eublepharidae began diversifying in the Cretaceous. Eublepharis split from its sister genera in Africa in the Palaeocene-Eocene, and began diversifying in the Oligocene-Miocene, with late Miocene-Pliocene cladogenesis giving rise to extant species. The current species diversity within this group is moderately underestimated. Our species delimitation suggests 10 species with four potentially unnamed divergent lineages in Iran, India and Pakistan. All 30 individuals of E. macularius that we sampled from the pet trade, which include diverse morphotypes, come from a few shallow E. macularius clades, confirming that lab and pet trade strains are part of a single taxon. One of the wild-caught haplotypes of E. macularius, from near Karachi, Pakistan, is identical to (10) pet-trade samples and all other captive populations are closely related to wild-caught animals from central/southern Pakistan (0.1-0.5 % minimum pairwise uncorrected ND2 sequence divergence).

Population Genomics of Wall Lizards Reflects the Dynamic History of the Mediterranean Basin.

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an "Out-of-Italy" expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

Checklist of Pseudoscorpions (Arachnida, Pseudoscorpiones) of Albania.

A checklist is presented comprising 42 species in 21 genera and 10 families of the arachnid order Pseudoscorpiones from Albania. Two taxa are identified only to the genus level -Neobisium Chamberlin, 1930 and Roncus L. Koch, 1873. The checklist was constructed according to available published data, material deposited in the Natural History Museum in Vienna and newly collected material. Individual families contained the following numbers of species: Chthoniidae (8), Neobisiidae (13), Garypinidae (1), Geogarypidae (2), Olpiidae (1), Cheiridiidae (1), Atemnidae (1), Cheliferidae (7), Chernetidae (7) and Withiidae (1). Sixteen species and four families (Garypinidae, Cheiridiidae, Atemnidae, Withiidae) are recorded for the first time in Albania -Chthonius jonicus Beier, 1931, C. rhodochelatus Hadzi, 1933,C. tenuis L. Koch, 1873, Ephippiochthonius serbicus (Hadzi, 1937), E. tuberculatus (Hadzi, 1937), Amblyolpium dollfusi Simon, 1898, Geogarypus italicus Gardini, Galli and Zinni, 2017, Apocheiridium ferum (Simon, 1879), Atemnus syriacus (Beier, 1955), Beierochelifer peloponnesiacus (Beier, 1929), Hysterochelifer cyprius (Beier, 1929), H. meridianus (L. Koch, 1873), Rhacochelifer peculiaris (L. Koch, 1873), Dinocheirus panzeri (C.L. Koch, 1837), Lamprochernes chyzeri (Tömösváry, 1883) and Withius piger (Simon, 1878). Moreover, A. dollfusi and G. italicus are newly recorded from the Balkans. Data on European distribution, habitat preferences and taxonomic considerations are reported. The Albanian pseudoscorpion fauna is compared with that of other Balkan countries.

Morphological and molecular data on tadpoles of the westernmost Himalayan spiny frog Allopaa hazarensis (Dubois & Khan, 1979).

Little is known about the life history, ecology, and distribution of the genus Allopaa (Dicroglossidae) and far less recent data are available about the larvae of this taxon. Here, we provide data on the larval stage of Allopaa hazarensis (Dubois & Khan, 1979) from northern Pakistan based on the examination of three tadpoles. Specimens were obtained from two sites in Buner, Khyber Pakhtunkhwa province, Pakistan. Morphological and genetic analysis (mtDNA and nDNA) confirmed the identity of the tadpoles as A. hazarensis. Tadpole characterizations were illustrated by detailed imagery. Basic measurements and details on oral apparatus provide relevant taxonomic characteristics to distinguish the tadpoles of this species from other spiny frogs. The illustration and description of the tadpole of A. hazarensis should facilitate the identification of this species in the field.

Relict groups of spiny frogs indicate Late Paleogene-Early Neogene trans-Tibet dispersal of thermophile faunal elements.

BACKGROUND: The Himalaya-Tibet orogen (HTO) presents an outstanding geologically active formation that contributed to, and fostered, modern Asian biodiversity. However, our concepts of the historical biogeography of its biota are far from conclusive, as are uplift scenarios for the different parts of the HTO. Here, we revisited our previously published data set of the tribe Paini extending it with sequence data from the most western Himalayan spiny frogs Allopaa and Chrysopaa and using them as an indirect indicator for the potential paleoecological development of Tibet. METHODS: We obtained sequence data of two mitochondrial loci (16S rRNA, COI) and one nuclear marker (Rag1) from Allopaa samples from Kashmir Himalaya as well as Chrysopaa sequence data from the Hindu Kush available from GenBank to complement our previous data set. A Maximum likelihood and dated Bayesian gene tree were generated based on the concatenated data set. To resolve the inconsistent placement of Allopaa, we performed different topology tests. RESULTS: Consistent with previous results, the Southeast Asian genus Quasipaa is sister to all other spiny frogs. The results further reveal a basal placement of Chrysopaa relative to Allopaa and Nanorana with an estimated age of ca. 26 Mya. Based on the topology tests, the phylogenetic position of Allopaa as a sister clade to Chaparana seems to be most likely, resulting in a paraphyletic genus Nanorana and a separation from the latter clade around 20 Mya, although a basal position of Allopaa to the genus Nanorana cannot be entirely excluded. Both, the placements of Chrysopaa and Allopaa support the presence of basal Paini lineages in the far northwestern part of the HTO, which is diametrically opposite end of the HTO with respect to the ancestral area of spiny frogs in Southeast Asia. These striking distributional patterns can be most parsimoniously explained by trans-Tibet dispersal during the late Oligocene (subtropical Chrysopaa) respectively early Miocene (warm temperate Allopaa). Within spiny frogs, only members of the monophyletic Nanorana+Paa clade are adapted to the colder temperate climates, indicating that high-altitude environments did not dominate in the HTO before ca. 15 Mya. Our results are consistent with fossil records suggesting that large parts of Tibet were characterized by subtropical to warm temperate climates at least until the early Miocene. They contradict prevalent geological models of a highly uplifted late Paleogene proto-Plateau.

Out of the blue: The first record of the genus Heremites Gray, 1845 (Squamata, Scincidae) from Pakistan.

The genus Heremites Gray, 1845 is endemic to the Western Palearctic region, containing morphologically similar species with a not well resolved taxonomy. The genus has a broad distribution from North Africa to Central Asia, with the only known record from northeastern Afghanistan. Three species are currently recognized in the genus with one, H. septemtaeniatus (Reuss, 1834), representing populations at the eastern edge of the genus range. During extensive fieldwork, we discovered H. septemtaeniatus from northwestern Pakistan and provisionally suggest that this population could be morphologically defined as H. septemtaeniatus transcaucasicus (Chernov, 1926). This important contribution to the knowledge regarding the family Scincidae in Pakistan, however, needs further investigation using an integrative approach.

Molecular and Morphological Data Confirmed First Record of Abbreviata kazakhstanica Markov and Paraskiv, 1956 (Spirurida: Physalopteridea) in Iran.

BACKGROUND: The genus Abbreviata (Spirurida: Physalopteridea) currently contains 47 species. Physalopteridae nematodes infect a large number of vertebrates, including mammals, birds, reptiles and amphibians. The current study is a report of the first morphological and molecular identification of A. kazakhstanica (Spirurida: Physalopteridea) in Pseudopus apodus in Iran. METHODS: Eleven road-killed P. apodus, were collected from, Iran during 2016-2018. The nematodes were isolated from stomach. After morphological study, the genomic DNA of the parasites was extracted using CTAB method. The DNA was used for PCR amplification of cytochrome c oxidase subunit I (cox1). The PCR products were sequenced, the sequence data were analyzed and multiple alignments were conducted using the Clustal Omega. RESULTS: After detailed microscopic examination, the A. kazakhstanica was identified. The cox1 sequences confirmed the species of helminth. The new sequences of A. kazakhstanica were submitted to GenBank under the accession number MK578751-2. CONCLUSION: Regarding the limited data on parasitological status of Iranian reptiles, more specific and comprehensive investigations are needed to identify the parasitic fauna.